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Abstract 

In this paper we study different algorithms for backward stochastic differential equations (BSDE 
in short) basing on random walk framework for 1-dimensional Brownian motion. Implicit and 
explicit schemes for both BSDE and reflected BSDE are introduced. Then we prove the convergence 
of different algorithms and present simulation results for different types of BSDEs. 
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1 Introduction 

Non-linear backward stochastic differential equations (BSDEs in short) were firstly intro- 
duced by Pardoux and Peng ([20J, 1990), who proved the existence and uniqueness of the 
adapted solution, under smooth square integrability assumptions on the coefficient and the 
terminal condition, and when the coefficient g(t,u,y, z) is Lipschitz in (y,z) uniformly in 
(t,u). From then on, the theory of backward stochastic differential equations (BSDE) has 
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been widely and rapidly developed. And many problems in mathematical finance can be 
treated as BSDEs. The natural connection between BSDE and partial differential equations 
(PDE) of parabolic and elliptic types is also important applications. It is known that only 
a limited number of BSDEs can be solved explicitly. To develop numerical methods and 
numerical algorithms is very helpful, both theoretically and practically. 

The solution of a BSDE is a couple of progressive measurable processes (Y,Z), which 
satisfies 



where B is a Brownian motion. Here £ is terminal condition and g is a generator. From 
[2D] , we know that when £ is a square integrable random variable, and g satisfies Lipschitz 
condition and some integrability condition, BSDE (pQ) admits the unique solution. 

The calculation and simulation of BSDEs is essentially different from those of SDEs (see 
[T5]). When g is linear in y and z, we may solve the solution of BSDE by considering its 
dual equation, which is a forward SDE. However for nonlinear case of g, we can not find the 
solution explicitly. Here we describe a software package that compute our numerical solutions 
for BSDEs with a convenient user-machine interface^. This package computes solutions of 
BSDEs, reflected BSDEs with one or two barriers as well as BSDEs with constraints. One for 
significant advantage of this package is that users have a very convenient interface. Any users 
who know the ABC of BSDE can use this package very easily. The input-output interface 
was also carefully designed. 

This paper is organized as follows. In Section 2, we introduce the discretization of 
BSDEs, then present implicit and explicit schemes for numerical calculation and consider 
their convergence. In Section 3, we show some numerical simulations. In Section 4, we 
consider reflected BSDEs with one barrier which is an Ito processe, by implicit reflected 
scheme, explicit reflected scheme, penalized explicit-implicit scheme and penalized explicit 
scheme, then we prove the convergence of these schemes. In Section 5, we apply penalized 
schemes to BSDEs with constraint on z and BSDE with solution y reflecting on a function 
of z. 

We should point out that there have been many recent different algorithms for computing 
solutions of BSDEs and the related results in numerical analysis, for example [2J, [3J, [I], [5], 
0, 0, 0, HQ, H3, H81, PU, [25], [2E], W\- In contrast to these results, our method uses 
very simple method. 

2 Numerical Schemes for Standard BSDEs 

Let (fl, JF, P) be a complete probability space, (B t ) t > be a 1-dimensional Brownian motion 
defined on a fixed interval [0,T]. We denote by {J~t}o<t<T the natural filtration generated 
by the Brownian motion B, i.e., Tt = cr{B s ; < s < t} augmented with all P-null sets of T . 

1 The study of simulations of BSDE has been started since 1996 in Shandong University, Mathematical 
Finance Laboratory directed by PENG Shige. First simulation was done by ZHOU Haibin, then following 
his works XU Mingyu worked on this software package since her master program(from 2000). This paper is 
a summary of almost all algorithms that have been used in the package. The algorithms for reflected BSDE 
with two barriers will be discussed in details in another paper. 




(1) 
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We consider for a fixed n G N, 

[*/*] T 

B n t := v^Ve", for all < t < T, 5 = -, 
f 1 J n 

where {e™}" =1 is a {1, — l}-valued i.i.d. sequence with P{s^ = 1} = P{£] = —1} = 0.5, i.e., 
a Bernoulli sequence. We set := cr{e", • • • , £™} and tj = Sj. 

Let g : [0, T] x R x R — ► R be a Lipschitz function in (y, z) uniformly of t, i.e., g satisfies 
for a fixed \i > 

yi, *i) -g(t,y 2 ,z 2 )\ < n(\y x - y 2 \ + \z x - z 2 \) (2) 
Vt G [0,T],\/(y 1 ,z 1 ),(y 2 ,z 2 ) G R x R. 

And <?(■, 0, 0) is square integrable. 

We will approximate a pair of real- valued (jF t )-progressively measurable processes (Y, Z) 
defined on [0,T] such that £[sup < 4 < T |y t | 2 ] + £[/ |Z t | 2 dt] < oo, which satisfies 

-dY t = g(t,Y t , Z t )dt - Z t dB t (3) 

with given terminal condition Yt = £ G L 2 (jF T ), where L 2 (JF T ) is the space of Tt measurable 
random variable satisfying E |£| 2 < oo. It is clear that Y has continuous paths. An existence 
and uniqueness theorem for equation ([3]) was established in [20], when the generator g satisfies 
<^ and g(-, 0, 0) is a square integrable. In many situations we are also interested in BSDEs 
of the following form: 

-dY t = g(t, Y t , Z t )dt + dA t - Z t dB t , t G [0, T], (4) 

where (A t )te[o,T] is an (.^-predictable RCLL process with almost surely bounded variation 
such that A Q = and £'[sup 0<t<T \A t \ 2 ] < oo. By the standard existence and uniqueness 
theorem for solutions of BSDE, for each given A and Y? = £ G L 2 (jF T ), there exists a unique 
pair (Y, Z) for equation (jlj). Here K has RCLL paths. We call the triple (Y,Z,A) a g- 
supersolution (resp. g-subsolution), if A is an increasing process (resp. decreasing process). 
It is called a (^-solution if A = 0. It is easy to check that, if both (Y, Z, A) and (Y, Z, A) are 
g-supersolutions on [0, T], then {Z, A) = [Z, A). Thus we often call Y a g-super(sub)solution 
(or ^-solution when A = 0) without specifying the related (Z, A). 

2.1 Implicit and Explicit Schemes for BSDEs 

We first give an assumption for discrete terminal condition £ n . 

Assumption 2.1 Consider £ which is Tt -measurable and £ n which is Q™ -measurable, such 
that 

Ei\Z\ 2 ] + sup E[\C\ 2 } < oo 

n 

and 

linage -T| 2 ] = 0. 
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Example 1 Set £ = $((S t )o<t<r), where $ : D[o,t] '-^ R one? satisfies Lipschitz condition. 
By Donsker's theorem and Skorokhod representation theorem, there exists a probability space, 
such that sup 0<t<T \B™ — B t \ — ► 0, as n — > oo, m L 2 (jF T ) ; since zs m L 2+<s . 5o £ n := 
$((i?")o<t<T); with!;, satisfies Assumption \2. 1\ 

The numerical solution of ([3]) is obtained by (Y t n , Zf) = (yj,Zj), t E + 1)5), 

5n = T. (w™, Zj)a<j< n is the solution of discrete BSDE which starts from = £ n . Our 
discrete BSDE on the small interval is 

y? = vj+i + ate, vh Z D 5 - z 7 £ l+^- ( 5 ) 

Then for given we want to find ^"-measurable (y™, 2"). The feasibility of this scheme 
for small S is due to the following easy lemma. 

Lemma 1 Lety™ +1 be a given Q™ +l -measurable random variable. Then, when 5 < 1/fi, there 
exists a unique £?" -measurable pair (y^, zf) satisfying (Tjp. 

Proof. We set Y + = y" +1 | E «_ i=1 and K_ = ?/5+i|e? +1 =-i- Both V + and Y_ are ^-measurable. 
Equation (jSJ) is then equivalent to the following algebraic equation: 

y] = Y + + g(t j ,y?,z?)6-z?V6, 

This is equivalent to 

z l = ^rS Y + ~ Y -) = -^E[yU^\Q]\- (6) 



2Vs ~ Vs 



and 



y] - gfavlzftS = \{Y + + K_) = E[y] +X \G]\. (7) 

Because g is assumed to be Lipschitz, the mapping Q(y) = y — g(tj, y, z™)8 is strictly mono- 
tonic: when SfM < 1, 

(e(y) - e(y'), y - 1/') > (l - fy) |y - yf > o. 

So there exists a unique value satisfying (j7j). □ 

This lemma shows a way to solve (jSJ), and we named this algorithm as 'implicit scheme'. 
In many cases, @ _1 cannot be solved explicitly. Thus we introduce the following explicit 
scheme by using E[y^ +l \Q^} to approximate y™ in g of © . We set Y£ = y™ = £ n and, starting 
from j = n — 1, solve in following reverse order, 



Then we get, 



V6 L ^+W^ J 2^ 
This explicit scheme is useful when g is not linear in y, for example <?(£, y, z) = sin(y). 
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Example 2 In pricing option, if deposit interest r and loan interest R are different, we get 
g(t, V, z) = ry + o-Qz + (R - r)(y - z)~ . 

Remark 1 To find g-super (sub) solution with an increasing process A as in Q), we need to 
consider the discretization of A, setting Aq = 0, A n - := J^i=o E[A ti+1 — A t .\Q™\. Since A is 
an increasing process, A™ is also increasing. Then instead of (TJP, we get 

y] = y- + i + g{tj, yj, z])8 + (A] +1 - A]) - z?$ +1 Vs, 

where Aj +l — A™ is Q -measurable. Then from implicit scheme we get 
y] = G-\E[y? +1 \g?} + (A] +1 -A])), 

«■ = ^m^UG?} 

And from explicit scheme, we get 



T s = m+il^+a^E^g^s + iA^-A]) 

^n _ 1 m -n c n , r m _ $+1 l e "+i =1 ~ ^ k+l^ 1 



In t/izs paper, we will not make special efforts to study the convergence of discrete g-super ( sub jsolution. 
Indeed, if we set 

y? = V? + A?,z? = %, 0<j<n, 

then (y n , z n ) is discrete solution of discrete BSDE with coefficient g(t, y, z) = g(t, y — A t , z). 
When A n — > A in certain sense, then we can get the convergence of (y n ,z n ) by (y n ,z n ), 
which is discrete solution of a classical BSDE. 

However in many cases, the increasing process A is not given, it is associated with (Y, Z) 
in order to keep (Y, Z) to satisfying certain condition, like reflected BSDE and constraint 
BSDE. We will discuss them later in this paper. 



2.2 Convergence Results for Numerical Schemes for BSDEs 

We set 

yn _ n yn _ n -yn _ -n 7" _ —n f\ < + < T 

1 t — y[t/s]^ Zj t — z [t/s]i r t — y[t/<5]> — Z[t/5], u ^ r ^ i , 

where (y™, <2™)o<j<n and (l/jN ^j)o<j<n are discrete solutions of ([5]) by implicit and explicit 
schemes, respectively. 

By Donsker's theorem and Skorokhod representation theorem, there exists a probability 
space, such that sup 0<t<T \B™ — B t \ — > 0, as n — > oo, in L 2 (jF T ), since is in L 2+<5 . Here 
L 2+s is the space of random variable <f> satisfying E[((f>) 2+S ] < +00. Then we have 

Theorem 2 We suppose that assumptions \2. 1\ hold and that g is Lipschitz in y and z. Then 
the discrete solutions {(Y n , Z n )}™ =l under the implicit scheme and {(Y n , Z n )}™ =1 under the 
explicit scheme converge to the solution (Y, Z) of (T3j] in the following senses: as n — > oo ; 

E[ sup \Y t n -Y t \ 2 + / \Z:-Z s \ 2 ds] ^0, (9) 

0<i<T Jo 
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and 

rT 

E[ sup \Y t n -Y t \ 2 + / \Z?-Z a \ ds]^>0. (10) 

0<t<T Jo 

The convergence (Q for this implicit scheme was obtained in 2001 by a profound result 
of Briand, Delyon and J. Memin [5], which can also be found in [6]. From these results, the 
convergence (llOp can be derived. Before proving fllOp . we first present following lemmas. 

Lemma 3 Let a, b and a be positive constants, 5b < 1 and a sequence (vj)j=i n of positive 
numbers such that, for every j 



v j 



+ a < a + b5 v^. 



i=i 

Then 

sup Vj + a < ae bT . 

This is a type of Gronwall lemma for discrete cases. The proof can be found in [Tj5], so 
we omit it. 

Lemma 4 We assume that 5 is small enough such that (1 + 2/x + 2fi 2 )5 < 1. Then 

i j= o 
where C c , g = (1 + 5fi)E[\Cf] + S 2 (^ 0, 0)5. 
Proof. From explicit scheme 

y] = y] + i + gfa, E[y] +1 \gj],z])5 - z]Vde j+1 . 

We have 

\y-\ 2 -\y] + i\ 2 = -\^\ 2 S + 2[y^g( tj ,E[y] +1 \g-),z]))5 (12) 

\-7rin \nn~\ -r-n\\2 ^2 



-2jfi%V8e j+1 + 2j%g(t j ,E$ +1 \g?),zf)8V5e j 



3+1 



Taking expectation and the sum for j — i, ■ ■ ■ , n — 1 yields 

n-l 

E\y?\ 2 < E\C\ 2 -J2 E \^\ 2S 



3 =i 

n-l 



3 =i 
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Since the last term is dominated by 

n— 1 

<^£Wi 2 (l + /. + 2 yU 2 ) + |^,,0,0)| 2 + / i|E[^ +1 |^]| 2 + -|^| 2 } 

re— 1 1 

< 5£ £{|^| 2 (1 + 2/x + 2^ 2 ) + |^(t if 0, 0)| 2 + -|^| 2 } + /i^le 



we thus have 

re— 1 re— 1 

2 



iwi 2 +^j>i^r<$ ^ x)b(t i ,o I o)i a tf+(i+/K5)Eiri 

re-1 

+5(l + 2 / i + 2/i 2 )£E|^| 2 

j=i 

Then by Lemma [3j we obtain 

snpE|yn 2 + -£^|^| 2 5<C r , 9 e( 



re-l 

3 (l+2 At +2/i 2 )T 



2 

j=0 



For (fTTj) . we recall (|T2l) . and take the sum for j = i, ■ ■ ■ , n — 1 and sup over j, then take ex- 
pectation. Notice that X)*- =0 ^^Voej+i an d J2)=o 9^j> E[V] + i\Gj], SVSe j+1 are both 
martingales with respect to Qf, we apply Burkholder-Davis-Gundy inequality for them with 
similar techniques as before, then get 



£[sup|$f] < cC^n + C^EW] 



11—1 
3=0 

< cC c gn +0^ sup E\y]\ 2 



With previous results, we obtain fill I) . □ 

Proof of Theorem [2j, The convergence of (Y n ,Z n ) to {Y, Z) is proved in [5]. To prove 
(ITU]) , the result for (Y , Z ), it suffices to prove as n — > oo, 

£[ sup |y t " - F/f + / T |Z? - Z n s \ 2 ds\ -> 0. (13) 

0<t<T Jo 

From (J5]) and (JHJ), we have 

Itf-lffl 2 = \y? +1 -y? +1 \ 2 -E\z?-z?\ 2 S (14) 
-\g{ti,v?,$) -g^m+xm^)?* 2 -2(y?-lff)W-%)V5e j+1 
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Then we take expectation and the sum over i from j ton- 1. With £ n - £ = 0, we get 

n-l 



i=j 



+2 £ J5[(y» - fl) {g{t v y» ^ ) - i^f??], 2?))]* 

1 71—1 71—1 

< -2^Ek n -^i 2 ]+v^Ebr- : 

i=j i=j 
n-l 

+2f,5Ej2\y? - y?\ ■ \y? - m+MW- 



V7\ 2 } 



Since Iff — E[y^ +l \Q"j] = g(tj, E\yj +1 \Qf\,Zj)5, the last term is dominated by 

71—1 71—1 

5J2(^ + my- - y?\ 2 + J2v 2 mtj,Eiyr; +1 \gj],r;)\^. 

i=j i=j 

But with (ITT]), the second term is bounded by C5 2 . We thus have 

r n—1 n—1 

E \v? -y]\ 2 + 2 E E W ~ ^ < (i + + 2/i 2 )5[^ em -v?\ 2 } + cs 2 

i=j i=j 

By Lemma [3], we get 

sup£ \y n 3 -y]\ 2 < C 5 2 e^ +2 ^ T . 

Then we reconsider square of the difference between the discrete solutions of implicit scheme 
and explicit scheme shown in ()14p . This time we first take the sum and sup,-, then take 
expectation. Using Burkholder-Davis-Gundy inequality and similar techniques, we get 



71—1 71—1 



3 i=j i= i 



z?-z?\ 2 } 



< CfjT sup E\y™ — y™ 



|2 



with previous results, ffTOj) follows. □ 



We now prove a more general result which will be useful in proving convergence results 
for schemes of reflected BSDEs. Consider the following BSDE 

-dY t = [gi(t,Y t , Z t ) + g 2 (t,Y t , Z t )]dt - Z t dB t , (15) 
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Here g\ and g 2 are both Lipschitz functions. Then we have the following implicit -explicit 
scheme to only replace by E[y^ +1 \Q^] in g 2 , 

y] = y? +1 + gx(tj,y], %)8 + g 2 (t J , \jff +1 \Q?], %)6 - z]e] +1 V8, (16) 

or, equivalent ly, 

y n 3 = E[yj + i\G?}+ 91(^,^,^)8 + g 2 {t h E[y? +1 \G?],z>})8, 

We also set Y t n = Z™ = zf t , s ^, < t < T. Meanwhile we consider the fully implicit 

scheme 

y n i = y] + x + gifo, yj, zj)5 + g 2 (t j , yj, zj)5 - zJe] +1 V~8, 
and let Y t n = y$ /S] , Z? = zf t/S] , < t < 1. 

Proposition 5 Under same assumptions of Theorem [H assume g\ and g 2 are Lipschitz 
functions. Let (Y, Z) be the solution of BSDE 173]) . Then as n — ► oo ; 

pT 

E[ sup \Y t n -Y t \ 2 + / \Z n s ~Z s ?ds\ (17) 

0<t<T Jo 

Moreover there exists a constant C 2 depending on T and fi 2 which is Lipschitz constant of 
g 2 , such that 

pT 

E[ sup \Y t n - Y t n \ 2 + / \Z n s - Z n s ? ds] < C 2 5 2 . 

0<t<T Jo 

The proof is similar to that of theorem [2] and we omit it. 

Remark 2 This scheme is very useful. For example, we will use it for penalization BSDE, 
which will be discusses in section 4-1- 



3 Simulation Results for BSDEs 

We consider the terminal condition Y? = £ which is a function of Bt'- Y? = £ = Q(Bt). In 
this case we set y,™ = £" = $(5^). It can be checked that our explicit schemes ([8]) (as well 
as the implicit scheme) will automatically derive 

j j 
y] := u(j, BJ S ) = u(j, $ = v(j, B%) = v(j, VsJ^e?). 

i=i i=i 

Since B™ s takes on j + 1 different values, the whole solution {?/", ^™}o<i<n-i is a 2- vector 
with wx (^ +1 ) values. For convenience, we set T = 1 in our simulation part. 

Applying the above numerical schemes, we have developed a Matlab toolbox for calculat- 
ing and simulating solutions of BSDEs. This toolbox starts with a Matlab figure window with 
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input area for generator g = g(t,y,z) and terminal function £ = where x stands for 

B T . Here g and $ can be any functions accepted by Matlab. These toolboxes can be down- 
loaded from |http:/ / 159.226.47.50:8080 /iam/xumingyu/English.jsp| by clicking 'Preprint' on 
the left side. 

Here we consider the case: g(t,y,z) = —5\y + z\, £ = = sin(|5i|). After in- 

putting these parameters of a BSDE, the numerical calculation for the BSDE are launched 
after clicking the button "calculate". When the toolbox indicates "the calculation is com- 
plete" , clicking any other button in button-area will produce different types of simulations, 
i.e., clicking "progress" will generate a figure displaying the dynamic evolution of backward 
calculation of states yj which starts from j = n and ends at j = 0. 

Clicking the button "B.M. and solution y" will produce the dynamic simulation of 
(t,B t ,Y t ), shown in Figure 1. 



center 




Figure 1: The solution surface with one trajectory 



Here a trajectory of Y t runs on a colored 3-dimensional surface represented u = u(t, x), where 
x stands for the space of Brownian motion B. 

Clicking "solution (y,z)" will generate another Matlab figure, displayed in Figure 2. This 
figure shows the 3-dimensional dynamic trajectories of (t, B t , Y t ) and (t, B t , Z t ) and, simul- 
taneously, 2-dimensional trajectories of (t, Y t ) and (t,Z t ). And there are two groups of 
trajectories on the figure 
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simulation of stochastic phenomena: y(t) 



simulation of stochastic phenomena: z(t) 




Figure 2: The trajectories of the solution 

We now compare some numerical solutions calculated by these algorithms: implicit 
scheme, explicit scheme and Monte-Carlo method in some particular situations. 

Case I. If g is a linear function (y, z): g{s,y, z) = by + cz + r. The solution Y of the 
BSDE is 

Y = exp((6 - \c 2 )T)E[ie^(cB T )} + ^[exp(6T) - 1]. 



Example 3 Set b = c = r = l, ^ = sind-Bpl). The numerical results obtained with the 
implicit and explicit schemes are shown in the following table: 



n 


100 


500 


1000 


2000 


5000 


~Vn 
J 


3.5106 


3.4916 


3.4879 


3.4866 


3.4859 


V'n 
J 


3.4171 


3.4716 


3.4785 


3. 4819 


3.4840 



The exact solution is expressed by Y = exp(|)_E[sin(|i?i|) exp(£?i)] +exp(l) — 1. We apply 
the Monte-Carlo method, with 10,000,000 samples, to calculate Yq. The result is Yq = 3.4850. 

Example 4 Set b = c = l,r = 0,^ = \Bt\- The numerical results obtained with the implicit 
and explicit schemes are: 



n 


100 


500 


1000 


2000 


5000 


Y 


3.1806 


3.1731 


3.1722 


3.1719 


3.1714 


Y 


3.0818 


3.1531 


3.1621 


3.1667 


3.1694 



Applying Monte-Carlo method with 10,000,000 samples to the exact solution 
Y = exp(|)E[|Si| exp(-Bi)], we get Y = 3.1710. 
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Case II. If g = ^z 2 , then we have the exact solution Y = In (£7 [exp (£)]). Since g does 
not depend on y, implicit schemes and explicit scheme give same results. 

Example 5 For £ = sin(|I?i|) ; applying the implicit scheme, we obtain: 



n 


100 


400 


800 


1000 


2000 


V'n 
r 


0.6249 


0.6253 


0.6254 


0.6254 


0.6255 



By Monte-Carlo method with 10,000,000 samples to the exact expressionY = ln(_E[exp(sin(|I?i|)]) ; 
we get Y = 0.6255. 



4 Reflected BSDEs 

4.1 Algorithms for reflected BSDEs with one barrier 



In this section, we discuss the algorithms for reflected BSDEs with one continuous lower bar- 
ier L. A solution o 

J Q T \Z S \ 2 ds + I-R'tI 2 ] < oo and 



rier L. A solution of such equation is a triple (Y, Z, K) on [0, T] satisfying E[sup 0<t<T \ Y t \ 2 + 



Y, = ( I g(s,Y s ,Z s )ds + K T -K t - [ Z s dB s 



t Jt 

T 



Y t > L t , < t < T, with f (Y t - L t )dK t = 0. 

Jo 

In [12], existence and uniqueness of the solution of this equation is proved when g satisfies 
Lipschitz condition (E2) and + /„ g 2 (t, 0, 0)dt + 

su Po<t<r(^) 2 ] < 00 • Here we consider 
the case when L t is an Ito process, i.e. L t = L + J* l s ds + f* a s dB s , < t < T and 
£ = $((-B s )o< s <t) satisfying requires of integrability, for convenience of discretization of 
processes. 

Remark 3 We call a progressively measurable process <p t is in space S 2 (0,T) ; if it satisfies 
£'[sup 0<f<T \4> t \ 2 } < oo. If a predictable process <p t is in space L^-(0,T), then it satisfies 

E[fy \<f) s \ 2 ds] < oo. And we define a space of Tt-measurable random variables which 

satisfies E[\^f] < oo, as lP{T t ), for (3 G R + . 

Following the same discretization introduced in section 2, we will approximate the so- 
lution of reflected BSDE. On the small interval [j5, (j + 1)5], the equation (|T8|) can be 
approximated by the discrete equation 

y] = vf+i + gfa&tys + q-tfefaVs, (19) 

> L],(y?-L])d] = 0, 

where d] = K t . +1 - K tj , and L] = L + Sjy^k + EtoWr+iV 7 ^ Here JTSJ) is called 
discrete reflected BSDE in [19], with terminal value £ n = $((Ei=o £ ?+i V / 5)o<j<n)- 
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Remark 4 When L t = ip(t,B t ) with ip £ C 1,2 ([0,T] x R), by ltd formula, we know that 
L t = L + / (J^ + \'^i)' l P{ s -, B s )ds + f* -§zip(s, B s )dB s . In fact, our algorithms are available 
for the case when the barrier L is a functional of Brownian motion, i.e. L t = ^/(t, (B s ) < s < t ), 
with its discrete version LVj., 5 i = ty(t[ t /s], (^1=0 £ k+i^")o<i<[t/s])- I n this section, we focus on 
ltd process in order to discuss the convergence of discrete solution. 

Suppose yj +1 is known, we try to find ^"-measurable (y™, z™, d") to satisfy f|T9|) . Set Y + = 
Vj+i\e* +1 =i and y_ = l«"+i=-i- From GB> we & et immediately = -±*E[y™ +1 e j+1 \G?) = 
-^(y+ — Y-)- Substitute it into the equation, our problem is changed to find (y™, d™) 
satisfying 

y] = E[yJ +1 \g^}+g(t v y^zJ)5 + d] } (20) 
Then we introduce two different schemes for this equation. 

Implicit reflected scheme. First, we present the implicit reflected scheme which is intro- 
duces by Memin, Peng and Xu in [19]. If we consider the mapping Q(y) := y — (g(tj, y, zj) — 
g(tj, LJ, z™))5, then for 6 small enough, we have 

(9(y) - e{y'),y - y') > (1 - 8fi \y - y'f > 0, 
i.e. Q(y) is strictly increasing with 0(L") = L™, so 

Q-\ y )>u;^y>q. 

It follows 

y] = e-\E[yJ +1 \g-}-g(t 3 ,L],z-)5 + d]), 
Notice that E[y] +1 \QJ\ = \{Y + + we get the results. 

Explicit reflected scheme Instead of solving the inverse of the mapping O, we replace 
by E | JFj 1 ] on the right side of (!20|) to get an approximal solution. Then it follows 

T S = m + i\^]+9(tj,E[y- +1 \gj]),z])5 + d r ;, (21) 
% = (m + t\^}+9(^m +1 \g-},z^)5-L]y. 

Substitute E^^Tj] = \(Y + + F_) into it, we get the results. 

Remark 5 Compare with the implicit reflected scheme, the explicit reflected scheme is much 
easier to compile programs for simulation. For example g(t,y, z) = sin(y). 



13 



Another important numerical method is via the penalization equations of reflected BSDE. 
In [12] , the authors introduced the penalization method to prove the existence of the solution. 
For p e N, the penalization equation with respect to the lower barrier L is 



Y? = £ + f g(s, YJ, Z*)ds + p J\y p - L s )-ds - f Z*dB s , (22) 



thanks to the comparison theorem for BSDE, we have Yf < Yf + , for p e N. Denote 
Kf = P J (Yg — L s )~ds. Then we know following results from [12]. 

Theorem 6 There exists a positive constant c independent on p, such that 

T 

\Z? - ZA 2 dt + sun \K P - KA 2 ] < ■ 
Whenp^ oo, we know Y p ^Y in S 2 (0,T) ; Z p -> Z in l£(0,T), K p -> K in S 2 (0,T). 



£[sup \Y t p -Y t \ 2 + I \Z P -Z t \ 2 dt + sup \K P -K t \ 2 ]< 

0<t<T Jo 0<t<T 



Numerical Penalization scheme By theorem [6j we know that the solution of reflected 
BSDE can be approximated by the solution of penalization equations ( 1221) . for some large p. 
Then on the small time [j5, (j + 1)5], we consider the following discrete penalized BSDE 



VT = V^i + 9it v y p -\ % n )5 + P(yf n - LfiS - z^4be 3 



If we have already known (yf+i, z j+i), then to solve (yf n , z™' p ) from above equation, we first 

get = ^Eiy^eUm = ^1 ~ where Y* = l^l^.Vf = ifflil^-x- 
Then y p ' n satisfies following equation 

VT = E[y^ x \g]\ + gfa yf\ + p(yf n - L])S. (23) 



There are two ways to find suitable y p,n '. One is penalization implicit scheme, i.e. to solve 
the equation: 

y P,n = (e^-\E[yf^\gj)) = {Q p )-\\{Yl + Y p )). 

Here 9 P is a mapping, Q p (y) = y — (g{tj, y, z p ' n ) + p(y — L™)~)5. Let d p,n = p(y p,n — Lu)~8. 

The other is implicit-explicit scheme, we only replace y p ' n of g in ([23"]) by Ely^'^J 7 ^]. 
Then we get, penalization explicit-implicit scheme, i.e. 



With E[y p '^ 1 \J : '"] = |(y?+i| e '? 1= i + Z/f+ile" i=-i)> resu hs follow easily. And we set d 
pfcj?" - L«)-5. 
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4.2 Convergence results of different schemes for Reflected BSDE 
with one barrier 

We first study the penalization scheme of reflected BSDE with one lower barrier. For pe- 
nalization implicit scheme, define Yf' n = y?tjg\, Z p ' n = zf^g, and Kf' n = Y^m%^rn ■ By 
Donsker's theorem and Skorokhod representation theorem, there exists a probability space, 
such that sup 0<t<T \B™ — B t \ — > 0, as n — > oo, in L 2 (jF T ), since Ek is in L 2+5 . For convergence 
of scheme, we have 

Proposition 7 Under assumptions \2.1\ and g satisfying Lipschitz condition. The sequence 
(Y p,n , Z p,n ) converges to (Y t , Z t ) in following sense 



r T 

lim lim E[ sup \Y t p ' n -Y t \ 2 + / \Z P ' n - Z s \ 2 ds] = 0, 

p->oon->oo <t<T Jo 



(24) 



and for < t < T, K p ' n — > K t in L 2 (jF t ), as n — > oo, p — > oo. 
Proof. Since 

£[ sup |F t p ' n - Y t \ 2 + / |Zf n - Z s \ 2 ds] < 2E[ sup \Y p ' n - Y t p \ 2 + / \Z p s ' n - Z s p | 2 (is] 

0<t<T Jo 0<t<T Jo 

+2£[sup \Y t p -Y t \ 2 + / \Z P -Z s \ 2 ds], 

0<t<T Jo 

by the convergence results of numerical solutions for BSDE and penalization method for 
reflected BSDE, Theorem |6l we know (1241) hold. For the increasing processes, we have 

E[(Kr - K t f] < 2E[(Kr - K p f] + 2E[{K P - K t f\. 

While for fixed p, 

K r = Yr-rr- [\(s,Yr,znd S + f z^dB:, 

Jo Jo 

K p = Y p -Y p ~ f g(s,Y p ,Z p )d S + f Z p dB s , 
Jo Jo 

from [B] Corollary 14, we know that f Z P ' n dB™ converges to J Q Z P dB s in S 2 (0, T), as n — > oo, 
then with Lipschitz condition of g and (12"4"1) . we get E[(K p,n — K p ) 2 ] — > 0, as n — ► oo. With 
convergence result of penalization methos, the result follows. □ 

Then we consider the penalization explicit-implicit scheme, note Y^' n = y^?g\, Z P ' n = 

and 'K P ' n = Y^m=o^m ■> ^ follows that 

Proposition 8 Under same assumptions of Proposition^ (Y t ' ,Z t ' ) converges to (Y t , Z t ) 
in following sense 

rT 

lim lim E[ sup \Y%' n -Y t \ 2 + / \Z P s ' n - Z s \ 2 ds] = 0, 

p-^oon->oo Q<t<T Jo 

with K^l' 71 —> K t in L 2 (jF t ) ; for < t < T , as n —>■ oo, p — > oo. 
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Proof. The convergence of (y t ' , Z f ' ) is a direct result of Proposition [5] and (I24I) . We 
consider the increasing process, notice that for < t < T, 

Jo JO 

with iff ' n = Y p ' n - Y t p ' n - j'l g(s, Y s p > n , ZP' n )ds + J* Z^ n dB n s , thanks to Lipschitz condition 
of g and the convergence of (Y P ' n , Z P ' n ), we get E[(Kf' n — 'K P ' n ) 2 } — > 0, as n — > oo, p — ► oo. 
With convergence results of penalization method, results follow. □ 

Now we study the convergence of reflected schemes. First for the implicit reflected scheme, 
denote Y t n = yf t/s] , Z? = zf t/s]l K? = Y}j=o d ^ for < t < T, from the results in pj, we 
know 

Theorem 9 (Theorem 3.2 in [19] ) Under assumption \2.1\ and (Hp for g, as n — > +oo, 

E[ sup \Y t n - Y t \ 2 ] + E I \Z? - Z t \ 2 dt -> 0. 

0<t<T Jo 

For the increasing process, we have 
Proposition 10 Fort e [0,7], E[(K t - K?) 2 } -> 0, as n -> oo. 
Proof. For t e [0, T], we have 

E[{K t - K?) 2 } < 3E[(K t - Kf) 2 ] + 3E[(K? - K p t ' n ) 2 } + 3E[(Kf' n - K?) 2 ] 
where K p is from penalization equation ([6]), and Kf' n is discrete solution of fl23|) . with 

R p,n = ^[t/5] d p,n_ g imilar ag Lemma 2 g . n we haye | F P,n _ y „|2] + £ j-T ^ _ 

Z™\ 2 dt < e ^ ,L , where C^ g L only depends on £ n , g, L and /i. Since 

*r = Yr-rr- 1 9 (s,Yr,zr)ds+ t z^b:, 

Jo Jo 

K n = Y n -Y t n - [ t g(s,Y s n ,Z:)ds+ tz:dB:, 
Jo Jo 

C R 

with Lipschitz condition of g, we deduce that E[(Kf' n — K^) 2 ] < . It follows 

s[(if t - K r t i ) 2 } < (c^ L + ci, L )^ + 3e[(k? - Krn 

Since Kf' n —>■ Kf in L 2 (jF t ) as n — > oo, for fixed p, we can choose n large enough to get right 
side very small. Then result of K n follows. □ 

Then we consider the convergence of the reflected explicit scheme. We set 

_ _ [t/5] _ 

j=0 

First as Lemma HI we have similar estimation of |L of reflected BSDE, given by (j2ip . 
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Lemma 11 We assume that 8 is small enough such that (2 + 2li + 6li 2 )5 < 1. Then 

where Cfn gL only depends on ll, E[\Cf], ES^fe. . ^ and £[sup,,-((L™) + ) 2 ] . 
Proof. Recall that for j = 0, 1, ...n — 1, (y™, z") satisfies 

y] = ^ +1 + ^,(i?[^ +1 |^]),^)5 + rf;-^ +1 v / 5, (25) 
Apply similar techniques of Lemma H] to fl25l) , we have 

+E\g(t j ,E^ +1 \S^)^W-E[(^) 9 \. 

In view of (y] - Lfjdj = and d- > 0, it follows 

E\y™\ 2 + E\z™\ 2 5 < E\y] +1 \ 2 + 2E[yJ +1 -g(t jl E[y] +1 \gj}^ 

+E[\g(t,,E[y] +1 \g-],zJ)\ 2 5 2 ] 

< E\y] +1 \ 2 + (S + 35 2 )E[\g{t^ 0, 0)| 2 ] + {U + ^ 2 5 2 )E[{z n 3 ) 2 } 

+5(1 + 2/x + V + 3// 2 5)£|^ +1 | 2 + 2£[(L J ") + • 3?] 
Notice that 3fi 2 5 < |, since 6/i 2 <5 < 1. Taking the sum for j = i, • • • , n — 1, it yields 



n—l 



3=1 



n—l „ n—l 

6 ~ ' ^"^'^ 1-2 

j+ll 

71-1 



< E \e? + (5 + 35 2 )E £>(f if 0, 0) | 2 ] + 6($ + 2/x + V)£ £ |gr 

n—l 

+aE[sup((L«) + ) 2 ] + - J E;[(^rf;) 2 ], 



where a is a constant to be decided later. Since d.- =Vj — Vj+i ~ 9{tj, {E\yj +1 \Q'j]) , Zj)S + 
zp] +1 V6, we get 

n—l n—l n—l 

E^i = Vi - ^ - TjSitj, (E[y] +1 \gj}),z^)6 + E^ £ i + iA 

j=i j=i j=i 

taking square and expectation on both sides, it follows 

n—l n—l n—l 

£[(E^) 2 ] ^ ^IVil 2 + 4 -£|6 n | 2 + W8TE Y^[\d(tj, 0, 0)| 2 ] + 12li 2 S E\fJ +1 \ 2 (26) 

n-1 



+45(3^ 2 + l)^|^| 2 . 
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Set a = 32, notice that 5(3/i 2 + 1) < |, then ^M±I2 < ^ we get 



7 E\y?\ 2 < 1 ^E\C\ 2 + & + 35 2 )Ej2[\9(t 3 ,0,0)\ 2 ]+32E[snp((L]) 
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n-l 



+ 5(| + 2^ + f / , 2 )E^|^ +1 | 2 . 



Then apply Lemma [3], in view of assumption that implies <5(| + 2/x + y/^ 2 ) < ^^ we obtain 

I — n I 2 



BupJ5[|j£f<q^. 



It follows from the estimations of Zj and that 

n—1 n—1 
j=0 j=0 

As Lemma HI using Burkholder-Davis-Gundy inequality and similar techniques, we get the 
results. □ 

Then we have following convergence result for explicit reflected scheme. 

Theorem 12 Under the same assumptions of Theorem^ the discrete solutions {{Y n , Z n )}™ =1 
of the explicit reflected scheme converges to the solution (Y, Z) of < f7%)) in the following senses: 
as n — > oo 

pT 

E[ sup \Y t n - Y t \ 2 + / \Z r s l - Z s \ 2 ds] -» 0. (27) 

0<t<T Jo 

Moreover E[sup 0<t<T (K t — ~K^) 2 } —> 0. 

Proof. Thanks to convergence results of Theorem [91 it suffices to prove 

E[ sup |y» - ?/f + / T |Z; - Z:\ 2 ds] -> 0. (28) 

0<t<T Jo 

Recall the implicit reflected scheme and explicit reflected scheme: for < j < n — 1, 

y] = y] +1 + g(t 3 ,E[y] +1 \gj]),zJ)S + d;-z]e] +1 V~S, 
Consider the difference, we have 

E \y?-y]f 

= E \y» +1 - ^Vi| 2 " ^ |^ - zf + 26E[(tf - y^Mtj, y» z«) - g(t 3 , E[y- +1 \g-],z]))} 
+2E[(y] - y])(d] - T 3 )} - 5 2 E[{g{t v y» z]) - g{t v E[y^ +1 \gj},z])f} 

-25E[(d] - S^fo, yl zj) - g( tj , E[y] +1 \gj),z]))) - E 



d--T 3 



2 
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in view of —a 2 — 2ab — b 2 = — (a + b) 2 < and 

~(y] - L])d] - (y] - L»)(d») 

< o. 

We take sum over j from i to n — 1, with £ n — £™ = 0, then get 

ra—l n— 1 

j=i j=i 

Now we are in the same situation as in the proof of Theorem [2j By similar methods, with 
Lemma [TTl and 

2n5E[\y^-y-\-\y--E[y] +1 \g-]\] 

= 2^E[\y] - f}\ ■ \y] - Ti + g{tj, E[y] +1 \G?\), $)6 + 3J 

< (2/i+ l)5E[\y]-y]\ 2 } +2^8E[8 2 \g(t 3 , E[y] +1 \g-}), $)\' + 



we obtain 

r n—1 n— 1 

where L is a constant only depends on £ n , 0, 0), /i and L. By Lemma [31 we get 

supE\y?-y]\ 2 <C5e^ +2 ^ T . 

j<n 

From ([H, it follows ([28]), which implies lim^oo 5 YTjZl E \zj - z™| 2 = 0. Then ([27} follows 
by using Burkholder-Davis-Gundy inequality, similar techniques and estimations results from 
Lemma [TTl In fact, we get 



n— 1 n—1 
|2 



z i I 



E[sup\y]-y]\ 2 ]<C ll 6j2E\y--y]\+Sj2 E 

j<n . . . . 

For the convergence of K™ , for < t < T, since 

[t/S] [t/S] 

K? = Y n -Y t n -J29(t J ,yJ,zJ)S + J2^] + i^, 
i=o i=o 

_ _ _ [t/6] [t/6] 

i=0 j=0 
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with Lipschitz condition of g and BDG inequality, we get 



[t/s] 




[t/s] 

+45(2/i 2 + c 2 )^|^-^ 



j=0 

[t/s] 




+2^j2\y?-y]\ 2 + sc c, 9 ,L- 

3=0 

From (1281) and convergence of to K t in L 2 (^), we obtain the convergence of to K t . 
□ 

4.3 Simulations of Reflected BSDEs with one lower barrier 

For calculation convenience, we consider the case when T = 1, and begin with y^ = £™ 
backwardly solve (yj, z 1 -, rf"), for j = n — 1, • ■ ■ 1, 0. Consider the amount of total calculation 
for most general case, we only treat a very simple situation: £ = (j)(Bi), L t = ip(t,B(t)), 
where <ft an d t/> are real regular functions defined on R and [0, 1] x R respectively. As for 
BSDE, we have also developed a Matlab toolbox for calculating and simulating solutions of 
reflected BSDEs. This toolbox is similar to the one for BSDE and can be downloaded from 
http://159.226. 47. 50:8080/iam/xumingyu/EnglishJsp| by clicking 'Preprint' on the left side. 

Here we consider following case: g(t,y,z) = — \y + z\, £ = $(-E>i) = 2sin(5 x ), L t = 
tf(t, B t ) = sia{B t + f ) - 2 and n = 400. 

After inputting the parameters, we run the calculation program using reflected explicit 
scheme, then get all prossible results of y. We may notice that at t = 1, £ > L\ does not 
always hold. But the numerical scheme still works as well. In fact, in such case the increasing 
process K as well as y has a jump of size (L\ — £) + at t = 1, which pushes the solution y t -, 
i.e. y n -i in our case, to stay above the barrier L. Then both K and y act as the terminal 
condition is (£ — L\) + + L\, which is always bigger than L%. 

Now we will see some properties of the trajectory of solution y in the Figure 3. In the 
upper portion of Figure 3, the below surface shows the barrier L in 3-dimensional, as well 
the upper one is for the solution y. Then we use programs to generate two trajectories 
of the discrete Brownian motion {B™' l )o<j< n , for % = 1,2, which are drawn on horizontal 
plane. The value of y^ii = 1, 2) with respect to these Brownian samples, are showed on the 
solution surface, and we use the fine vertical line to give correspondence between two group 
of trajectories of y and B. The remainder of the figure shows respectively the trajectories 
of the force = J2k=o^k' 1 (* = 1? 2) corresponding to the value of y^ l {% = 1,2), and 

n.i t n,i i ■ i o\ 
fj ~ L 3 ( Z = 1>2 )- 

In the upper portion we can see that there is an area where two surfaces (the solution 
surface and the barrier surface) stick together. When the trajectory of solution y™ goes into 
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n.i 



this area, the force K- ' will push upward. Indeed, if we don't have the barrier here, y 
intends becoming smaller than the reflecting barrier L™' 1 , so to keep y™' 1 being no less than 
L™' 1 , the action of forces K™ ,% are necessary. Comparing these two trajectories, we can see 
that one trajectory, noted as K™' 1 pushes upwards the corresponding trajectory of solution 



y" 1 ' , while the other one noted as Kj' 2 , 
the trajectory y™' 2 with respect to K™' 2 does not 



keeps zero, since y j ' goes into the sticking but 



j 



center 




0.1 



0.05 



I 



0.2 



0.4 0.6 
t 



0.8 1 0.2 

Figure 3: The solution on surface 



0.4 0.6 

t 



0.8 



Compare the two sub-figures, which are below the main sub-figure, we can easily find 

L™'* takes the value 0; but the converse is not 



out that the Kj' 1 only increases when y™' 1 

n,i rn,i g j^n4 



always true, when y™' 1 — L T - ,% = 0, KJ' 1 does not necessary increase. 

About this point, we can also see in Figure 4. This figure shows two groups of 3-dimensional 

i i2 ) and, simultaneously, two 



Z J )(*= 



dynamic trajectories (tj, B™'\ l^ n ' l )(i=i,2) and {tj,B r - 
groups of 2-dimensional trajectories of (tj, Yj 1 ' 1 ) ( i=1)2 ) and (tj, Z™ 



)(i=l,2). 



For remainder sub- 



figures, the above-right one is for the trajectories K- ' 



1,2), and while the below-left 



one is for y 



n,t 
j 



1,2), then comparing the these two sub-figures, as in Figure 3, we 



can see clearly the relation between K\ 



1,2) and y n { 1 - L) 



1,2). Moreover one 



trajectory, noted as K n,x as well as 7/ n,1 jumps at t 
the barrier. 



n,i 

1, since its terminal value is less than 
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trajectory of y(t) in 3-D trajectory of z(t) in 3-D 

simulation of stochastic phenomena: Ki t I 




Figure 4: Simulation of y, z, K 



Then we list out some numerical results for the reflected scheme and explicit implicit 
penalization scheme, and we can see that as the penalized parameter p converge to infinity, 
t/q' 71 converge to y$. Consider the same parameters as above: f(y,z) = —\y + z\, £ = 
$(-Bi) = 2sin(i?i), L t = ^(t,B t ) = sm(B t + |) - 2. Then as the following tablet showing: 

n = 400, reflected explicit scheme: y% = —0.6430, 



penalization scheme: 

n = 1000, reflected explicit scheme: 
penalization scheme: 

n = 2000, reflected explicit scheme: 
penalization scheme: 
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Remark 6 For BSDE with two reflecting barriers, we introduced also reflected implicit and 
explicit scheme as well as penalization schemes. The proofs of convergence and simulations 
results can be found in [2^]. 
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5 T-constrained BSDEs 

In this section, we consider a g-supersolution (Y, Z, A) with constraint (Y t , Z t ) G Y t of the 
following form: 

Y t = £+ I g{s,Y s , Z s )ds + A T - A t - I Z s dB s (30) 



with dr(Y t , Z t ) = 0, a.e., a.s., where T is a nonempty closed subset of M x I and dr is the 
distance function of T, i.e., dr(y,z) = iaft y i )Z r^ e p{\y — y'\ + \z — z'\}. It is clear that T is a 
Lipschitz function 

\d r (y, z) - d v (y' } z')\ < \y - y'\ + \z- z'\. 

Such a g-supersolution (Y, Z, A) is called a T-constrained g-supersolution. 

As before we assume that g satisfies Lipschitz condition (j2J) and that assumption 12.11 
holds for £. From [2T], we have the existence of smallest solutions for (|3"0"|) : 

Theorem 13 If there exists at least one T- constrained g- supers olution of Ultyl , then the 
equation admits a smallest T- constrained g- supers olution (Y,Z,A). Moreover, (Y, Z) is the 
limit of the following sequence of penalization solutions: 

Yt P = i+ f g{sX s P ,Zl)ds + p f d T {Y! , Zl)ds - f Z?dB s , (31) 
Jt Jt Jt 

in the sense of 

lim E [ T [\Y t - Y t p \ 2 + \Z t - Zff}dt = 0, 1 < (3 < 2. 

This smallest T-constrained (7-supersolution is called a gr-solution. And such equation 
can be considered as a BSDE with singular coefficient g? := g(t, y, z) + 00 • d r (y, z). It easy 
to check that when £ + e L°°(jF T ), and there exists a large enough constant C such that for 

y>c 

g(t,y,0)<C o + n\y\, and (y,0)eT, 

there exists a T-constrained ^-supersolution of (l30j) (see Peng and Xu [23]). Then by Theorem 
[T5| a (^r-solution exists. In this section we will work under these assumptions. We now derive 
a numerical scheme applying convergence results in Theorem [TSJ 



5.1 Constraint on Z 

First we consider the case when constraint is only on process Z and invariant in t , i.e. 
T is a close subset in R. And we require Z G T, i.e. dr(Z t ) = 0, a.e. a.s.. After same 
discretization for BSDEs introduced in Section 2, for each positive number p we have the 
following penalization discrete equation on small interval [j8, (j + 1)5] 

VT = 1/J& + 9(tj, V? n , zT^ + Vdv{z^)5 - zf n V~6e j+1 , (32) 
with discrete terminal condition: := £ n . 
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Now we need to find a way to find ^"-measurable (y p,n , z™' p ) to satisfy (1521) with [y p ^x, z ^+i)- 
It is easy to get zf n = -^E[y p ^ x e] +1 \J^] = ^(yffi| £ y=i - yffil^=-i)- Substitute it into 
( |32l) . it follows a equation of y p ' n as 

So apply the implicit scheme for BSDE in Section 2, we get 

y p ' n = Q-\E[y p ^\g]]+pd T {z p ^), 
where Q(y) = y — g(tj,y, z p ' n )5. While the explicit method gives 

The interesting point here is that the penalization of z p,n with respect to z is not directly on 
z p ' n , it act on y n,p to influence z p,n . 
We have 

Theorem 14 (Convergence Theorem) Define 

vP, n — „,p> n 7P, n — v p ' n \ 7P ' ri — r-P> n ~7 P)U — -^P> n 

I t — y[t/s\-> Zj t — z [t/5\-> 1 1 — y[t/5\-> — z \t/s\- 

Here y p,n and y P ' n , < j < n, can come from either implicit scheme or explicit scheme 
respectively. Under assumption ^. 11 and g satisfying Lipschitz condition. Then 

lim lim E[ / \Y p ' n - Y s \ 2 ds + / \Z P ' n - Z s f ds] = 0, 1 < j3 < 2. 

Proof. By Theorem [131 f° r an Y £ > 0, there exists po > such that for each p > p , 



E[[ \Y P -Y s \ 2 ds 2 + [ \Z P S -Z s fds]<e. 
Jo Jo 



Moreover, by Theorem [2J for implicit scheme, we have as n — > oo 



r-T 

E[ sup \Y P0 ' n - Y po | 2 + / \Z po ' n - Z P0 \ 2 ds] -> 0. 

0<t<T Jo 

For explicit scheme, the result follows from 

sup E[\Yf' n - Y p "\ 2 + / \Z P ; ' n - Z po \ 2 ds] -> 0. 

0<t<T Jo 

□ 

To illustrate calculation and simulation in our software package, we consider the case 
T = [a, b] with a < < b. Then dr(z) = (z — a) - + (z — b) + . The default setting is 
g(t,y,z) = -2 |y + z\ - 1, £ = with a = -0.5, b = 0.8, p = 20 and n = 400. The 
surface u = u p,n (t, x) and v = v n ' p (t, x) are given with dynamic simulation Y p,n = u p ' n (t, B™) 
and Zf ,n = v p,n (t, B™) as shown as the upper part of Figure 5. The lower part of the figure 
displays the simultaneous action of the process A P ' n = pEj<[t/^r( 2 f n )- The effect of 
increases in the A P ' n when Z p,n is less than —0.5 and larger than 0.8 are clearly shown in 
Figure 6. But it seems that the solution is still too sensitive to the choice of p and n. If 
p\fS > 1, then the numerical solution will explode. 
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5.2 BSDE reflected on process Z 

Now we consider another special case, when the constraint is dr(y,z) — (y — <p(z))~, in 
other words, we require y — <p{z) > 0. After the same discretization of the time interval, we 
have following discrete penalization equation for some p large enough, on the small interval 
[j5,(j + l)5],0<j<n-l 

vT = vf+i + 9te, + p( W) - vT) +5 - # n 8e S+ i. 

Similarly, zf n = ^[y^e^Tf] = ^(yf^\ £ n =1 - ^; n 1 | e n = _ 1 ). Then y p ' n satisfies 

Set Q(y) = y-(9(^ y ,z^+p(^ n )-y) + S), with E[y p ^] = |(#AI^+lfAI^=-x), 
then our implicit scheme is given by solving following equation 

yf n = e- 1 {E[ffi 1 \g»)). 

Meanwhile, we have also explicit-implicit scheme, which is 

As previous subsection, we have convergence results of these two schemes. 



Theorem 15 Define Y t p ' n = y p ^ 6] , Zf' n = zj^ and Ff™ = y p ;J s] , Z P '" = z p [J s] . TIk-ii we 
have, for 1 < (3 < 2, 



lim lim E[ \Y p > n -Y s \ 2 ds + \Z p s ' n - Z s f ds] = 0, 

p^oo n-»oo J Q J Q 
pT pT 

lim lim E[ \Y P s ' n -Y/ds+ 1^' - Z s f ds] = 0. 



p— >oo n— >oo 

Proof. The results follow from Theorem [TBI and Proposition [5j so we omit the proof. □ 

Now we do simulations by explicit-implicit scheme. We consider the case g = —2 \y + z \ — 
1, £ = <p(z) = 1.25 x z, with penalization parameter p = 10, and discretization number 
n = 400. In figure 7, we see the surface of solution Y p ' n with a trajectory of Y p ' n on 
the surface in upper portion, while in two lower subfigures there presents the trajectory 
AP,n _ p^j < , t ,J K y P j ,n — (f)(z p,n ))~S and yf n — 4>{z p ' n ) corresponding to the one on the surface. 
We can see that when y p,n — (p(z P,n ) is positive, the penalization term will not work to the 
process y*j' n - About this point we can see more clear in Figure 8, which presents trajectories 
of y p ' n , z p ' n , A p ' n and y P,n — 4>(z p,n ) in 3 or 2-dimensional subfigures. 
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Figure 7: The solution surface of penalization BSDE with dr(y, z) = (y — 4>(z)) 




0.5 1 0.5 1 0.5 1 



t t t 

Figure 8: A trajectory of solutions of penalization BSDE with d r (y, z) — (y — <fi(z))~. 
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